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ABSTRACT 

Autonomous Underwater Vehicles (AUVs) are self-propelled robotic platforms that carry a suite of 
instruments and sensors which enable the vehicle to navigate the depths of the ocean using its onboard power. Due to 
the limited power onboard the vehicle, efficient streamlining is an important design consideration for maximum 
endurance of a cruising type AUV. Torpedo-shaped hull forms are very common among cruising type AUVs, which 
consist of a nose section, a cylindrical mid-body section, and a tail section. This paper discusses different torpedo-shaped 
bare-hulls formed by changing the nose and tail sections to identify which combination of nose and tail renders the most 
hydrodynamically efficient hull form. A torpedo-shaped AUV bar e-hull for which experimental results were available is 
selected as the base model and CFD analysis is performed and validated. With the validated CFD model , the different 
shaped hull configuration is changed by modifying its nose and tail profiles. A total of seven configurations were 
analyzed for three different speeds. The least drag shaped configuration is identified and recommended for further 
inclusion of control surfaces. 
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Nomenclature 

F D Drag force (N) 

X Length (m) 

R Maximum radius of nose profile (m) 

N Nose index 

0 Semi-tail angle 

U Speed (m/s) 

C Tail length (m) 

T Time (s) 

0 Total control volume bounded 

A Total length of nose (m) 

V Volume of displacement of the hull (m 3 ) 

C dv Volumetric Drag Coefficient 

1. INTRODUCTION 

An Autonomous Underwater Vehicle (AUV) is a self-propelled robotic vehicle controlled and piloted by 
an onboard computer, which is capable of maneuvering in 3D space without the supervision of an operator. Under 
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most environmental conditions the vehicle can follow precise pre-programmed trajectories. Sensors onboard the AUV 
sample the water column, providing the ability to make both spatial and time series measurements. Data collected by an 
AUV is automatically geospatially and temporally referenced and normally of superior quality. Multiple vehicle surveys 
can increase productivity and insure adequate temporal and spatial sampling thereby provide a means of investigating the 
coherence of the ocean in time and space [1]. AUVs find application in various fields such as research, oil exploration, 
defense purpose, etc. Underwater vehicles fall into two basic categories: (a) Manned and (b) Unmanned. The manned 
versions include submarines and passenger carrying submersibles. The costs and risks associated with the manned 
underwater vehicles are more that in the recent past decades there has been an obvious trend toward unmanned vehicles. 
Unmanned underwater vehicles include towed vehicles, ROVs and AUVs. A towed vehicle normally does not have any 
propulsion system and is generally connected to vessel by a cable. An ROV is a tethered vehicle but has thrusters for 
maneuvering in 3D space and generally controlled by an operator from a surface ship. In addition, they are powered 
through tethers and transfer of data is established from the vessel operator to receive commands during its operations. 
As far as AUVs are concerned, they are not connected by tethers or cables for its operation during a mission. The fact that 
AUVs are untethered brings in the necessity for onboard power on an AUV. Hence, the endurance of an AUV is limited by 
the battery power they carry onboard. An autonomous underwater vehicle is typically designed to operate at a nominal 
speed range of 0.4 to 2 m/s [2]. Since, the availability of battery power on board is limited, to minimize the powering 
needs, a streamlined hull-form is generally preferred for a cruising type AUV. This reduces the resistance experienced by 
the hull in the fluid medium. There are different hull forms but the most common hull-forms preferred for cruising type 
AUVs are: i) Torpedo shaped profile ii) Low drag profile, Figures 1 and 2 

A low-drag profile has the advantage of a more streamlined flow around the body and hence the drag is less. 
However, constraints like fixed size (lack of modular flexibility) and costlier fabrication constitute the difficulties [12]. 

„ I 


I-Tail section- 1 I-*—Variable payload-1 i-Nosesection'-j 

Figure 1: Torpedo Shaped Profile (Nose, Cylindrical Mid-Body, Tail), [7] 



Figure 2: Low Drag Streamlined Profile (Radius Varies along the Length), Sea-glider [11] 

A torpedo-shaped AUV as shown in Figure 1 on the other hand, is modular in design consisting of a forward nose 
section, a mid-body and an aft tail section. The easy to fabricate cylindrical mid-section and the feasibility of modular 
interchange outweigh the demerits such as increased drag and structural weight. As hull geometry largely influences the 
ability of an AUV to efficiently maneuver in the water, the focus of this study is to compare the effect of hull geometry 
(profile shape) on hydrodynamic efficiency. In this paper, a torpedo-shaped hull is selected and numerical analysis has 
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been carried out using CFD techniques by modifying the bare hull with a different nose and tail configurations to obtain 
the hydro dynamically efficient (least drag) profile. 

2. METHODOLOGY 


The bare hull selected is the standard profile proposed by D. F. Myring with a cylindrical midbody [9]. The 
bare-hull geometry presented in [6]is taken as the base model, see Figure 3. The equations defining the nose and tail 
profiles are given below, 


r(x) 



(t 5 ) 2 ]" <*«> 


r(x) 


d 3d 
2 2c 2 


tanO 


0 - hf + [J- 0 - Iff, (Tail) 


( 1 ) 

( 2 ) 


Where, r = maximum radius of nose profile, a = total length of nose, x= length, n- nose index, a= semi-tail 
angle, c= tail length etc. Variables such as nose index ( n ), semi-tale angle (0°), were identified responsible for the fullness 
of both nose and tail sections respectively. This model had a semi-elliptical nose profile which confirmed to a shape 
described by Eq.(l) when the Myring parameter n takes a value equal to 2 and the tail is termed as Huang tail [4]. The 
methodology adopted is to obtain the least drag profile through numerical investigation using CFD. The experimental 
results available [9] is taken as the reference data. 



Figure 3: Bare-Hull Geometry (all Dimensions in mm) 

Nose Shape Myring, (n=2) with Huang Tail 

In the present study, a total of seven different profiles consisting of four different noses and three tails has been 
configured and analyzed. These different configurations were obtained by changing the parameters n and 00 in the Myring 
equation. The configured profiles of nose and tail shapes are shown in Figures 4 and 5 respectively. 

n = 1.25 

n — 2 

n = 3 

Hem i-spherical 

Figure 4: Configured Nose Profiles using Eq (1) and Hemispherical Shape 
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The analysis has been carried out in three phases. In the first phase, the selected combination and the available 
experimental results from the literature are validated with CFD. In the second phase, with the validated CFD solver 
parameters, the nose shape of the same model is being replaced by a different combinations. The AUV hull models of 
different nose profiles with Huang tail [4] is shown in Figure 6. The low drag nose profile is found from these 
combinations, i. e., n = 3, and then the tail of the AUV models are replaced by combination of tail profiles as shown in 
Figure 7. Finally, a least drag hull profile is obtained with a combination of least drag nose and tail profile amongst all 
models. Total seven profiles are tested in both the stages at different speeds (0.6, 1.0 and 1.4 m/s) and drag coefficients of 
all the models are compared and tabulated. 


L D = 10 



Myring 25 deg tail 


Hemispherical nose 



Figure 6: AUV Models with Modified Nose Profiles Figure 7: AUV Models with Modified Tail Profiles 


3. NUMERICAL MODELING 


To investigate the single-phase flow of fluid around the bare hull geometry, RANSE solver has been used for the 
numerical investigation using the Siemens STARCCM+ module, considering the viscous effects. The equation for RANSE 
is given below: 

jf ipPq ,dv+ f r p(p(v — v b ).ndS = TVcp. ndS + J n q dlV dV (3) 

The parameters Q is total control volume bounded by control surface is given by \\t. ‘t’ is time, p is the density of 
the fluid medium, v is velocity vector of fluid. 

The computational domain around the bare hull geometry of AUV is taken as shown in Figure 8 where the 
velocity inlet is at 10D from the nose end in bow, pressure outlet is at 39D from tail end of the bare hull, wall boundaries 
are at 5D from both in top and bottom and along port and starboard sides. D is the maximum diameter of the bare hull i. e. 
0.14m. 
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Figure 8: Computational Domain 

The domain is divided into three regions for 3D grid generation, the region differs in term of the size of the cells 
and the associated growth ratio of cell dimensions perpendicular to the surface of the AUV hull. The mesh generated 
around the hull is shown in the Figures 9 and 10. The boundary conditions used for the computational domain are shown in 
Table 1. 



Figure 9: Volume Mesh Generation around the Hull 




Figure 10: Volume Mesh Generation on AUV Bare Hull 
Table 1: Boundary Conditions for the Computational Domain 


Parts 

Boundary Condition 

AUV bare hull 

Wall no-slip 

Inlet boundary 

Velocity inlet 

Outlet boundary 

Pressure outlet 

Exterior boundary 

Wall-slip 


4. RESULTS AND DISCUSSIONS 

The model [9] is first put into numerical investigation using CFD. The results at three speeds (0.6, 1, 1.4 m/s) are 
computed and compared with the experimental work for validation. The drag force is obtained from the CFD solver and 
the results are plotted for the coefficient of drag. The following empirical formulae are used for obtaining the drag 
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coefficient, 

c » = < 4 > 

Where C dv is the volumetric coefficient of drag, F D is the drag force, V is volume (m 3 ) of displacement of the hull 
and U is the speed (m/s) and p is the density (kg/m 3 ) of the fluid medium. Next, the hulls are modified with the nose and 
tail profiles keeping the L/D ratio as 10 for the uniformity and comparison. Due to the changes in the nose and tail profiles, 
there is minimal variation in overall volume among the AUV models. The CFD results obtained has only a maximum 
percentage error deviation of 1.45%. Henceforth the domain size and boundaries applied here is considered as the reference 
for further investigations with different configurations of AUV hulls. The validation of the CFD results and the 
experimental work from the literature are shown in Figure 11. Table 2 shows the percentage error deviation of the CFD 
results with the published literature. 



Figure 11: Comparison of CFD Results with the Available Experimental Results 


Table 2: Error Comparison of Experiment with CFD Results 


Velocity U (m/s) 

Experimental Results 
(Literature) 

CFD Results 

Deviation % 

c dv 

c dv 

0.6 

0.0489 

0.0481 

1.45 

1 

0.0419 

0.0421 

0.58 

1.4 

0.0389 

0.0387 

0.43 


4.1 Comparison Results for Different Nose Profiles 

With the validated computational domain, the drag estimation is carried out for different nose profiles with Huang 
tail combination for the selected speeds. The comparison of these results on the drag coefficient with the modification in 
the nose is shown in Figure 12 and tabulated in Table 3. 
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Figure 12: Comparison of the Drag Coefficient for Different Nose Geometry 


Table 3: Drag Coefficient Values for Different Nose Profiles 



Drag Coefficient Values for Different Nose Profiles (C dv ) 

Base Model (n=2) and 
Huang Tail 

n =1.25 

© 

ro 

II 

a 

Hemi-Spherical 

0.60 

0.0482 

0.0550 

0.0477 

0.0510 

1.00 

0.0421 

0.0430 

0.0417 

0.0440 

1.40 

0.0387 

0.0390 

0.0377 

0.0370 


It can be seen that for the nose with n=1.25 and Hemispherical shapes, the coefficient of drag is maximum at 
lower speeds and reduces when the speed increases. In the case of the hemispherical nose, the separation occurred in the 
forward part due to a sudden change in the geometric profile. But when the speed increases, it is observed that the drag 
reduced only at the maximum considered speed. From this analysis, it can be stated that the Myring nose, n=3 profile 
produces less drag constantly moving from lower to higher speeds. An error comparison has been made and tabulated in 
Table 4 with the base model and the selected nose profile n = 3 for further studies by incorporating the tail sections. The 
deviation(%) shows that the obtained drag coefficient results from the CFD results for the model with nose profile n=3 is 
better than all the other profiles. 

Table 4: Error Comparison of C dv of Base Profile (Literature) and n=3 


Speed 
(m/s) 

Volumetric Drag Coefficient Values (n) 

Deviation 

% 

C dv for n=2 

C dv for n=3 

0.60 

0.0482 

0.0477 

1.05 

1.00 

0.0421 

0.0417 

0.96 

1.40 

0.0387 

0.0377 

2.65 


4.2 Comparison Results for Different Tail Profiles 

The model with nose n=3 and Huang tail is considered as the least drag profile combination for the bare hull. It is 
further required to investigate these models by changing the tail profiles of the hull. The optimum tail profile is obtained by 
changing the value of semi-tail angle (0) in the Myring tail equation. Three profiles are considered. The models are then 
tested with the inclusion of tail profiles to identify the hydrodynamically low drag profile for the bare hull. The CFD 
simulation has been carried out with the modified hull with different tail profiles. The comparison of the results is shown in 
Figure 13 and the results are tabulated in Table 5. 
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Figure 13: Comparison of Drag Coefficient for Different Tail Profiles 


Table 5: Drag Coefficient Values for Different Tail Profiles 


Speed (m/s) 

Coefficients of Volumetric Drag (C dv ) 


For Different Semi-Tail Angle (0) 

Huang tail 

(0)=O # 

(0)=15° 

(0)=25° 

0.6 

0.0482 

0.0501 

0.0476 

0.0497 

1 

0.0421 

0.043 

0.0415 

0.0426 

1.4 

0.0387 

0.0392 

0.0382 

0.0391 



^ Hiiang tail 
Myring (0)=O 
Myring (0)=15 
Mating (0)=15 



0.6 


0.S 1 

Speed (m/s) 


1.2 


1.4 


It can be seen that the tail profile with 15° semi-tail angles is the optimum tail with Myring, nose n=3. The results 
from CFD analysis shows that, the combination of Myring nose (n=3) and Myring tail semi-tail angles (15°) is having the 
minimum drag. The comparison of the final AUV model with the modified tail, i. e., 15° semi-tail angles are tabulated and 
presented in Table 6. The percentage deviation indicates that tail profile with (0) = 15° is having lesser drag coefficient 
values than Huang tail. 


Table 6: Comparison of Least Drag tail with Huang Tail 


Speed (m/s) 

Volumetric Drag Coefficients (C dv ) 

Base Model (n=2) and Huang Tail 

(0) = 15° 

Deviation % 

0.6 

0.0482 

0.0476 

1.3020 

1.0 

0.0421 

0.0415 

1.3374 

1.4 

0.0387 

0.0382 

1.3419 


5. CONCLUSIONS 


In the present investigation, a numerical analysis using the CFD technique has been carried out to obtain the low 
drag profile for the selected AUV models. This is an initial research work for developing a biomimetic AUV model with 
higher capabilities of hovering and surveillance for military purpose. The numerical results are first compared and 
validated with the published research to obtain standard settings in the CFD tool. Various nose and tail profiles have been 
generated with the Myring equation and further the CFD analysis has been carried out. The Myring nose with n = 3 with 
semi tail angle of 0, 15° is found to be the lowest drag profile among the selected hull forms. In the future, the hull would 
be modified with biomimetic propulsion with the control surface and experiments to be carried out to validate the 
developed model. 
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